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Abstract 

We establish a general principle for the tomographic approach to quantum state recon- 
struction, till now based on a simple rotation transformation in the phase space, which allows 
us to consider other types of transformations. Then, we will present different realizations of 
the principle in specific examples. 



1 Introduction 

The tomography is well known in the field of medicine where it is extensively used for images 
reconstruction in diagnostic systems. It is based on the possibility of recording transmission 
profiles of the radiation which has penetrated a living body from various directions. A collection 
of these distributions, by means of computer assisted mathematical procedures |1| , allows to obtain 
the desired physical density distribution giving insight into the body. 
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An approach, based on this principle, to get the Wigner function of a quantum system was 
first proposed by J. Bertrand and P. Bertrand ||. 

In quantum optics one has the opportunity of measuring all possible quadratures of an e.m. 
field by means of the homodyne detection, so that the tomography can be easily implemented. In 
fact, recently, Vogel and Risken pointed out that the homodyne detected marginal distribution 
is just the Radon transform (or "tomography") of the Wigner function. By inverting the Radon 
tranform one can obtain the Wigner function and then recover the state; this is the basis of the 
method proposed by Smithey et al. [Q. 

The density matrix elements, in some representations, can also be obtained by avoiding the 
Wigner function and then the Radon transform ]5|, || [J]. In particular in Ref. [J], the system 
density operator was expressed as a convolution of the marginal distribution of the homodyne 
output and a kernel operator. 

Anyway, also at quantum optical level the tomographic reconstruction relies on the possibility 
of performing measurements of observables obtained by means of transformations belonging to the 
group 0(2). Here, starting from this transformation, we would establish a general principle which 
include other types of transformations, i.e. other schemes for the quantum state measurement. 
This will be done in Sec. 2, while in the following we shall consider different realizations of 
this principle. In particular, in Sec. 3 we present the symplectic tomography, in Sec. 4 the 
latter is compared with the homodyne tomography and in Sec. 5 we discuss the photon number 
tomography. 

2 The tomographic principle 

In the usual Optical Homodyne Tomography the observed quantities are the quadratures = 
q cos (J) + ft sin obtained as mixtures of position q and momentum p by means of a rotation g in 
the phase space 
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The quadrature histograms w(x, 0), also called marginal distributions, are projections (Radon 
transformations) of the Wigner function [§] 

w(x,(f)) = J W(q cos <p — p sin <p, q sin (p + p cos <p) dp . (2) 

From the set of histograms w(x,(p) the Wigner function itself (hence the quantum state) can be 
reconstructed via the inverse Radon transform as was shown in Ref. 0. 
On the other hand, the marginal distribution w(x,<p) results as 

w(x, 0) = (xtlplxt) = (q\G(g)pG' 1 (g)\q} , (3) 

where \x^) are eigenkets of quadrature operators and G(g) is the unitary group representation for 
the transformation g. In this case 

(4) 

Thus the above well known state reconstruction procedure could be generalized into the fol- 
lowing principle: given a density operator p and a group element g, one can create different types 
of tomography if, by knowing the matrix elements (x\G(g)pG~ 1 (g)\x) , from measurements, is able 
to invert the formula expressing the density operator in terms of the above "marginal" distribu- 
tion (the x may denote either continous or discrete eigenvalues). For the inversion procedure one 
can use the properties of summation or integration over group parameters. The only problem is 
mathematical one to make the inversion and/or physical one to realize the transformation G(g) in 
laboratory. 

It is worth noting that other tomographic-like approaches already known |], [Kj can be taken 
back to this principle, but we would now put our attention on two other particular realizations. 

We would also remark that the marginal distribution function w depends on one variable 
and is determined by some extra parameters. Since we are treating the latter in the same way 
of variables, in the follows, we often use the steatment marginal distribution instead of set of 
marginal distributions. 



G{g) = exp 



2 + 2 . 
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3 Symplectic tomography 
3.1 Formalism 

Let us consider the quadrature observable X as a generic linear form in position q and momentum 
V 

X = fiq + up + 5 (5) 

with /i, u, 5 real parameters (their physical meaning will be discussed later on), then it is possible 
to get the density matrix elements from the marginal distribution, avoiding the evaluation of 



the Wigner function as an intermediate step. For this pourpose we start from a well known [] I 
representation of the density operator 

P = f ^W(q,p)f(q,p) (6) 

where the Wigner function W(q, p) is a weight function for the expansion of the density operator in 
terms of the operator T(q, p), which is defined as the complex Fourier transform of the displacement 
operator D(£) £[T 



T (q.,,) = I ^ e ^-l)v-i^- P )u (7) 



Following the lines of Ref. we may write the marginal distribution w for the generic 

quadrature of Eq. (R) as 



w[x, fl, v) 



-^-^W(q,p)^^, (8) 



with x = X — 5. By means of the Fourier transform of the function w one can then obtain the 
relation 

W(q,p) = (2it) 2 z 2 w(z, -zq, -zp) (9) 

where — zq, —zp, z are the conjugate variable to /i, z/, x respectively and the Fourier transform w 
has the property 

w(z, -zq, -zp) = — w{l, -q, -p) . (10) 

z z 
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It is worth remarking that in this case the connection between the Wigner function and the 
marginal distribution is simply guaranteed by means of the Fourier transform instead of the 
Radon one. 

We now present a specific example of the measurable probability distribution to; to this end 
we consider the Schrodinger cat state of the type discussed in Ref. |J, i.e. 

| a + ib) + | a — ib) 



{2[1 + cos(2a6) exp(-26 2 )]}V2 ' 
with a and b arbitrary real numbers. Hence in this case we have O 



;n) 



w(x, fl, v) 



2x1/2 
IT J 

x < cosh 



1 



jJ? + V 2 



1/2 



[l + cos(2a6) exp(-26 2 )] 



exp 



(x — /ia) 2 + b 2 v 2 



Avb{x — /j,a) 
fl 2 + v 2 



cos 



2b{2yx - a(fi 2 - v 2 )) 
fi 2 + v 2 



(12) 



The result of Ref. can be reproduced as partial case of this formula. 

To get an invariant expression of the density operator in terms of the marginal distribution w, 
we insert Eq. (^]) into Eq. and expressing the marginal distribution in terms of the Fourier 
transform, we have: 

dqdp 



P 



or, in a compact form 



(27^ 



T(q,p) / dxdjjidv z w(x,[i,v)e 



P 



-izx+ifizq+iuzp 



dxd[idv w(x, fi, v)K^ v , 



(13) 



(14) 



where the kernel operator Kn, v is given by 



Kiiv = J_^ e -i*» e -^("-v)a t e ^("-NM)a e -4o* a +^) 



2tt 



(15) 



The fact that K^ v depends on the z variable as well (i.e. each Fourier component gives a selfcon- 
sistent kernel) shows the overcompleteness of information achievable by measuring the observable 
ofEq. @. 

If, and only if, the kernel operator is bounded every moment of the kernel is bounded for all 
possible distributions w(x, fi, v\ then, according to the central limit theorem, the matrix elements 



of Eq. ( fn|) can be sampled on a sufficiently large set of data. This is only possible in the number 



or coherent states basis, as can be evicted from the expressions fll5D , analogously to the results of 
Ref. 0. 



We named the developed procedure "Symplectic Tomography" [[yj, since in this case the 
"marginal" distribution is obtained by using a symplectic transformation g belonging to the sym- 
plectic group ISp(2, R) 

q \ / cos 4> sin 



Q 
P 
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P 
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sin 



cos 



A 

A- 1 



For this transformation one has 



H = A cos (j); v = A 1 sin 



0. 



(16) 



(17) 



Thus, for the realization of the scheme, the element g is the product of squeezing and rotation 
operators. It means that for our scheme the representation operator is 



G(g) = exp 



MM 



exp 



(qp + pq) 



A related tomographic approach based on symplectic group was formulated for short pulse mea- 
surements in Ref. |T^]. In the two dimensional case the symplectic group reduces to SL(2, R), but 
the formalism is still valid at higher dimensions for transformations belonging to Sp(2n, R) (n > 1) 



13 1 , and the scheme might be generalized to other Lie groups different from the symplectic one. 



3.2 Applications 

The quadrature of Eq. @ could be experimentally accessible by using for example the squeezed 
pre-amplification (pre-attenuation) of a field mode which is going to be measured (a similar method 
in different context was discussed in Ref. |0|). In fact, let a be the signal field mode to be detected, 
when it passes through a squeezer it becomes a s = a cosh s — aJ sinh s, where s is the squeezing 
parameter |16|]. Then, if we subsequently detect the field by using the balanced homodyne scheme, 
we get an output signal proportional to the average of the following quadrature 



(19) 
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where (p is the local oscillator phase. Eq. ( |T9| ) can be rewritten as 

E{<j>) = -j= (a[e -< * cosh s - sinh s] + a f [e^ cosh s - e -< * sinh s\) , (20) 

which coincides with Eq. (|J) if one recognizes the independent parameters 

|U = [cosh s — sinh s] cos 0; z/ = [cosh s + sinh s] sin <fi . (21) 

The shift parameter 6 has not a real physical meaning, since it causes only a displacement of the 
distribution along the X line without changing its shape, as can be evicted from Eq. @. So, 
in a practical situation it can be omitted. To be more precise, the shift parameter does not play 
a real physical role in the measurement process, it has been introduced for formal completeness 
and it expresses the possibility to achieve the desired marginal distribution by performing the 
measurements in an ensemble of frames which are each other shifted; (a related method was early 
discussed in Ref. |17| ). In an electro-optical system this only means to have the freedom of using 
different photocurrent scales in which the zero is shifted by a known amount. An experimental 
method based on an observable similar to that of Eq. fl5|) was also proposed in Ref. ||18|| . 

Our approach can be extended to multimode systems and, in the two-mode case it is inter- 
esting to use the connection between the Wigner function and the marginal distribution of only 
one quadrature for the case of heterodyne detection |L3| (particulary used to detect multimode 
squeezed states). 

Beside that it is worth noting that the marginal distribution w(x,p,u) introduced in the 
symplectic tomography consents an alternative formulation of the quantum dynamics in terms 



of classical probability distributions |19[. In fact it has the characteristics of reality, positivity, 
normalizability and measurability. Then, we can derive the evolution equation for the marginal 
distribution function w using the invariant form of the connection between the marginal distri- 



bution and the density operator, given by the formula (14j). From the equation of motion for the 
density operator 

dtp=-i[H,p], (22) 
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we obtain the evolution equation for the marginal distribution in the form 

J dxdf/,dv [w(x, fi, u, t)K^ v + w(x, /i, v, t)/^} = (23) 

in which the known Hamiltonian determines the kernel through the commutator 

I^ v = i[H,K,^]. (24) 

The obtained integral-operator equation can be reduced to an integro-differential equation for the 
function w in some cases. To do this we represent the kernel operator in the normal order 
form (i.e. all the momentum operators on the left side and the position ones on the right side) 
containing the operator K^ v as follow 

: 4„ := K{p) : K^ u : V{q) (25) 

where TZ(p) and V(q) are finite or infinite operator polynomials (depending also on the parameters 
H and v) determined by the Hamiltonian. Then, we calculate the matrix elements of the operator 
equation (|23|) between the states (p\ and \q) obtaining 

J dxdfidv {w(x, /i, v, t) + w(x, fi, u, t)1l{p)V{q)} (p\ : K^ v : \q) = . (26) 

If we suppose to write 

n{p)V{q) = U(p, g) = E E c nA^ ^ v)v n q m , (27) 

n m 

due to the particular form of the kernel in Eq. (|I5"D, Eq. (|26| ) can be rewritten as 

J dxdfidv {w(x, ft, v, t) +w(x, fi,u, t)Tt(p,q)} {p\ : K^ v : \q) — 0, (28) 

where p, q are operators of the form 

fid a \ ( % d v \ 

+ H "={-lT» + 2 z )-< (29 > 

while z, in the space of variables x, /i, v should be intended as the derivative with respect to x, i.e. 

Z ~ ^ (30) 



and when it appears in the denominator is understood as an integral operator. Furthermore the 
right arrow over IT means that, with respect to the order of Eq. (P7|) , the operators p and q act 
on the right, i.e. on (p\ : K^v '■ \q). Under the hypothesis of regularity of w on the boundaries, 
we can perform integrations by parts in Eq. (p8|) disregarding the surface terms, to get 



J dxdfidiy {w(x, fi, u, t) + w(x, //, u, t)tl (p, g)| (p\ : K^ v : \q) = , (31) 

where now tl means that the operators p, q 

fid a \ fid v \ 

'=(~;*~H *={- z T»-ry (32 > 

act on the left, i.e. on the product of coefficients c n ^ m (— z, fi, v) with the marginal distribution w. 
Finally, using the completness property of the Fourier exponents given by (p| : K^^ : \q) we arrive 
at the following equation of motion for the marginal distribution function 

d t w + wtl (p, q) = 0. (33) 

Thus, the distribution function which depends on extra parameters obeys a classical equation 
which preserves the normalization condition of the distribution. In this sense we always can reduce 
the quantum behaviour of the system to the classical behaviour of a set of marginal distributions. 
Of course, this statement respects the uncertainty relation because the measurable marginal dis- 
tributions are the distributions for individual variables. That is the essential difference (despite of 
some similarities) of the introduced marginal distribution from the discussed quasi-distributions, 
including the real positive Q-function, which depends on the two variables of the phase space and 
is normalized with respect to these variables. This our approach to the quantum dynamics remind 
the formulation of quantum mechanics without the wave functions made in Ref. [|20|| . In reality, 
it is not guaranteed (from a rigorous mathematical point of view) that given the initial condition 
for the function w there exist a solution of Eq. (|33"1), but we demonstrate fl9| the existence of 



an invertible map connecting the Schrodinger equation with Eq. (|33|) , so there is a one to one 
correspondence between their solutions. Then, any solvable quantum systems should admit an 
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evolutionary solution for w. For example, a harmonic oscillator in an initial coherent state \aco), 
gives [fL9] 



w(x, p, l>, t) 



x exp 



1 



2 , 2 

A* + z/ z 



-1/2 



exp 



-2 1 ck 1 2 — —5- + 2a/2— (Re{a } cos t — Im{a } sin t) 



v 



—x + v^RejctoH/U sin t + v cos t) + v / 2Im{a }( zy sin t — p cos t) 



p 2 + u 2 \z/ 



(34) 



4 Comparision between symplectic tomography and ho- 
modyne tomography 

In this section we connect the discussed realization of the tomographic principle with the usual 
one, i.e. the homodyne tomography. The relations among symplectic tomography and optical 
homodyne tomography are also discussed by Wunsche [21]. To get the density operator in terms 
of the marginal distribution, analogously to Eq. (|l4f) , we can start from another operator identity 



such as [11 



J 7T 

which, by the change of variables \i = — valuta, v = v^2Rea, becomes 

p = — [ dfidu Tr{pe~ lji }e lji = — [ dpdv Tr{pe^ £ }e J 
2n J z7r J 



(35) 



(36) 



where x = X — 5. The trace can be now evaluated using the complete set of eigenvectors {\x)} of 
the operator x, obtaining 

Tr{pe~ i& } = f dx w(x,p,u)e~ ix (37) 



then, putting this one into Eq. (|3"6"P, we have a relation of the same form of Eq. { TA ) with the 
kernel given by 

1 ... 1 ..... 

(38) 



1 ... 1 ..... 
K u = e~ lx e lx = e -«BgW+wp 



2tt 2tt 

which is the same of Eq. (^) setting z — 1. It means that we now have only one particular Fourier 
component due to the particular change of variables (the most general should be zp = -\/2Ima 
and zv = VzReai). 
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In order to reconstruct the usual tomographic formula for the homodyne detection J7[ we need 
to pass in polar variables, i.e. /i = — r cos <fi, v = —r sin <fi, then 

x — ► — ritf, = — r[gcos (f> + psin 0]. (39) 

Furthermore, denoting with a;^ the eigenvalues of the quadrature operator x^, we have 

Ti{pe~ i£ } = Tr{pe ir& *} = J dx^ w^x^^e 1 ™* (40) 

and thus, from Eq. ( |36| ) 

p = Jd(pdx (j) w(x < j > ,(f))k^ (41) 

with 

K 4> = — I dr re*^*-*** (42) 
which is the same of Ref. [0. Substantially, the kernel of Eq. is given by the radial integral 



of the kernel of Eq. (RST), and this is due to the fact that we pass from a general transformation, 



with two free parameters, to a particular transformation (homodyne rotation) with only one free 
parameter, and then we need to integrate over the other one. 

5 Photon number tomography 

Here we investigate another possible realization of the tomographic principle. More precisely we 
derive an invariant relation connecting the density operator of the radiation field with the number 
probabilities |22" |. So that in this case the measured observable has a discrete spectrum. 



Indeed Wallentowitz and Vogel ||| proposed the s-parametrized Wigner function reconstruc- 
tion similar, in its essence, to the present one by using direct photon counting, and contemporarely 
an analogous scheme was adopted by Banaszek and Wodkiewicz |24| . 



We named the state retrieval by direct photon counting " Photon Number Tomography" since it 
is made by just detecting the number of photons at a given reference field and then scanning both 
its phase and its amplitude; differently from the usual homodyne tomography where a marginal 
distribution is recorded by homodyne measurements and then scanning only the phase. The 

11 



method can be used either when the output beam is mixed with a reference field at a beamsplitter 
as in Refs. I 



|, P4 |, or in a physical situation similar to the one proposed by Brune et al. [^5| 
for the generation and measurement of a Schrodinger cat state, allowing thus the possibility of 
reconstructing its density matrix and, more generally, the possibility to reconstruct cavity QED 
field state when there is not a direct access to the field |2l 



As claimed in Ref . [TlJ , the generalized version of Eq. 

d 2 a. 



is 



P 



( —W{a,s)f{a,-s) 

J TV 



(43) 



where the s-ordered wheight function W(a, s) may be identified with the quasiprobability dis- 
tributions Q(a), W(a) and P(a) when the ordering parameter s assumes the values —1,0,1 
respectively; while the operator T represents the complex Fourier transform of the s-ordered 
displacement operator l)(£,s) = -D(£)e s '^' 2 / 2 , which can also be written as |26[ 



T(a, s) 



—D(a) fS 



D-Ha). 



(44) 



1- s v ' Vs - 1 

On the other hand the weight function W^a, s) is the expectation value of the operator T(a, s 
TTJ , i.e. W(a, s) = Tr{pT(a, s)}, then we obtain 

s + l^" 



d 2 a ~ _2 

1 



7T " 1-S VS-1 



(n\D(a)pD- L (a)\n)T(-a, -s) 



(45) 



Thus, it becomes clear that the weight function W(— a, s) is related to the ability of measuring the 
quantity (n\D(a)f)D~ l (a)\n) by scanning the whole phase space [27[|, i.e. by just varying a. Now, 
we may consider one mode of the radiation, whose state p one wants to reconstruct, contained 
inside a cavity and, immediately before the photon number measurement, a coherent reference 
field is " added" |25|, |28j , so that we may recognize 



w{n, a 



Tr {D(a) pi)' 1 {a)\n)(n\} = (n\b(a)pi>-\a)\n) 



(46) 



as the probability to detect n photons after the injection of the reference field a. Again, as an 
example of the above marginal distribution we consider the state of Eq. (|TT|) , for which we have 

1 



w(n, a) 



2[1 + cos(2a6) exp(-6 2 )]n! 



-\a+a+ib\ 2 



\a + a + ib\ 2n + e 



-\a+a-ib\ 2 



\a + a — ib\ 



2n 
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+ 2e- |a+a+4b|2/2 - |a+a -* |2/2 Re{(a + a + ib) n {a* + a + ib) n } 



(47) 



The addition process we are considering, following Ref. p5l , "is quite different from the 



combination of fields produced by a beam splitter, which mixes together distinct modes coupled 
to its two ports and introduces vacuum noise even in the absence of any classical input field" . We 
are indeed describing a much simpler field amplitude superposition mechanism, discussed in the 
Glauber's pioneering work pSfl. The photon number distribution (fill) results as the projection 



of the field state p over a displaced number state f29fl . In this realization the transformation Q 
is represented by the displacement operator, but here the "marginal" distribution is a discrete 
probability. 



The photon counting could be made either by means of atoms , as in the case of microwave 
cavity field, or by direct detection of the outgoing optical total field. Furthermore, setting 

2 f s + 1\ n - 

K s (n,a) = - T(-a,-s), (48) 



1-s Vs- 1 
we have, from Eq. ([45]) 

p = Y] /— w(n,a)K s (n,a) . (49) 

n=0 J n 

Thus, analogously to Eq. (|i~4]), we may assert that a density matrix element can be experimentally 
sampled if there exist at least one value of the parameter s inside the range [—1,1] for which the 
corresponding matrix element of the kernel operator K s is bounded. From Eqs. fl48|) and ([44]) , one 
immediately recognizes that this is possible for s e (—1, 0], in the number, coherent and squeezed 
representations and not in the coordinate basis. For s = —1, in Eq. (HI) n = only survives, 



however as was shown in Ref. |TlJ, the operator T becomes singular and can be used to construct 
an arbitrary density matrix when it is only weighted with a well behavied function. It means that 
w{0,a) = Q{a) in Eq. ©. 

We would stress the fact that the Photon Number Tomography, in contrast to optical homodyne 
tomography, does not need of sophisticated computer processing of the experimental data, the 
quantity measured in the experiment is proportional to the quasiprobability distribution at a given 
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phase space point. In particular our scheme becomes especially useful in the intracavity optical 



tomography where other similar schemes [£3|, £J] are not applicable; in fact it can be adopted in 



a situation in which the photon number is measured indirectly using a sequence of atoms passing 
through the cavity [j30fl , with the quantum efficiency, determined only by the duration of the 



measurement process (i.e. the length of the sequence), that could be very high p5 |. In that case, 
the scheme has also the advantage of being QND. Furthermore, it results suitable for cavity QED 



characterization, like Ref . |31]] , allowing the reconstruction of nonclassical states as well, which are 
extremely sensitive to losses and then their detection seems prohibitive by means of an outgoing 
field as in Refs. J23], ^3 . 



The principle of Photon Number Tomography turned out to be useful also for the reconstruction 
of the motional quantum state of a trapped atom f32|| . 

Finally, a method for direct sampling of the density matrix in the Fock basis, based on the 
photon number tomography, is implemented in Ref. [p3|1 , resulting as an improvment since it has 
been shown that in this case it is sufficient to vary only the phase of the reference field, keeping 
its amplitude constant. 

6 Conclusions 

We have shown that different methods for the state reconstruction could be ascribed to the same 
principle. In that sense we have provided a sort of unification among different (analytical) ap- 
proaches. At the same time we have given the physical interpretation of the presented realizations 
and the new observables introduced offer the possibility of using various measurements, different 
from the usual ones, to determine the state of the system under study. The relations between 
various measurable probabilities are provided in a forthcoming paper [ 04j| . 
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